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Ultra-high-energy, > 10^^ eV, cosmic-ray and high energy, ~ 10^^ eV, neutrino pro- 
duction in GRBs is discussed in the light of recent GRB and cosmic-ray observations. 
■ Emphasis is put on model predictions that can be tested with operating and planned 
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1. Introduction 



The origin of GRBs, bursts of 0.1 MeV — 1 MeV photons lasting for a few seconds, 
Q^l remained unknown for over 20 years, primarily because GRBs were not detected prior to 
p • 1997 at wave-bands other than 7-rays |[l|]. The isotropic distribution of bursts over the sky 
^ ■ suggested that GRB sources lie at cosmological distances, and general phenomenological 
considerations were used to argue that the bursts are produced by the dissipation of the 
kinetic energy of a relativistic expanding fireball (see [2] for review). 
^ I Adopting the cosmological fireball hypothesis, it was shown that the physical conditions 
' in the fireball dissipation region allow Fermi acceleration of protons to energy > 10^°eV 
[^,^, and that the average rate at which energy is emitted as 7-rays by GRBs is comparable 
to the energy generation rate of UHECRs in a model where UHECRs are produced by a 
cosmological distribution of sources Based on these two facts, it was suggested that 
GRBs and UHECRs have a common origin (see [5] for a recent review). 

In the last two years, afterglows of GRBs have been discovered in X-ray, optical, and 
radio wave bands 0. Afterglow observations confirmed the cosmological origin of the 
bursts, through the redshift determination of several GRB host-galaxies, and confirmed 
standard model predictions of afterglows that result from the collision of an expanding 
fireball with its surrounding medium 0. These observations therefore provide strong 
support for the GRB model of UHECR production. 

In this review, UHECR and neutrino production in GRBs is discussed in the light of 
recent GRB and UHECR observations. The fireball model is briefly described in §2.1, and 
proton acceleration in GRB fireballs is discussed in §2.2. Implications of recent afterglow 
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observations to high energy particle production are discussed in §3. Model predictions 
are shown to be consistent with the observed UHECR spectrum in §4. Predictions of the 
GRB model for UHECR production, that can be tested with future UHECR experiments, 
are discussed in §5. High energy neutrino production in fireballs and its implications for 
future high energy neutrino detectors are discussed in §6. 

2. UHECR from GRB fireballs 

2.1. The fireball model 

In the fireball model of GRBs, a compact source, of linear scale tq ~ 10^ cm, produces 
a wind characterized by an average luminosity L ~ lO^^ergs"^ and mass loss rate M = 
L/rjc?. At small radius, the wind bulk Lorentz factor, F, grows linearly with radius, until 
most of the wind energy is converted to kinetic energy and F saturates at F ~ ~ 300. 
Variability of the source on time scale At ^ ro/c ~ 1 ms, resulting in fluctuations in 
the wind bulk Lorentz factor F on similar time scale, lead to internal shocks in the 
expanding fireball at a radius F^cAt. These shocks reconvert part of the kinetic 
energy to internal energy, which is then radiated as 7-rays by synchrotron emission of 
shock-accelerated electrons. The 7-ray flux is observed to vary on time scale t^ar ~ 



As the fireball expands, it drives a relativistic shock (blast wave) into the surrounding 
gas. At a radius r ~ F^cT, where T ~ 10 s is the wind duration, most of the fireball energy 
is transferred to the surrounding gas, and the flow approaches self-similar expansion. 
The shock driven into the ambient medium at this stage continuously heats new gas, 
and accelerates relativistic electrons that produce by synchrotron emission the delayed 
radiation, "afterglow," observed on time scales of days to months. As the shock-wave 
decelerates, the emission shifts with time to lower frequency. 

2.2. Fermi acceleration in GRBs 

The observed GRB and afterglow radiation is produced by synchrotron emission of 
shock accelerated electrons. In the region where electrons are accelerated, protons are 
also expected to be shock accelerated. This is similar to what is thought to occur in 
supernovae remnant shocks |0. Internal shocks are generally expected to be "mildly" 
relativistic in the fireball rest frame, i.e. characterized by Lorentz factor 7^ — 1 ~ 1, since 
adjacent shells within the wind are expected to expand with Lorentz factors which do 
not differ by more than an order of magnitude. We therefore expect results related to 
particle acceleration in sub-relativistic shocks to be valid for the present scenario. In 
particular, the predicted energy distribution of accelerated protons is dNp/dEp oc E~^. 

Two constraints must be satisfied by fireball wind parameters in order to allow proton 
acceleration to Ep > 10^^ eV in internal shocks [Q: 



in order for the proton acceleration time ta to be smaller than the wind expansion time, 
and 



Ti/T^c ~ At. 



(1) 




(2) 



3 



in order for the synchrotron energy loss time of the proton to be larger than ta- Here, 
r = 3001300, At = lOAtioms ms, Ep = lO^o^p^ao eV, = lO^^^^ ^^g^g/g jg ^^e 7-ray 
luminosity, is the fraction of the wind energy density which is carried by magnetic 
field, 47rr^cr^(i?^/87r) = C,bL, and C,e is the fraction of wind energy carried by shock 
accelerated electrons. 

Eqs. (|ip and (j^) imply that protons may be accelerated in a GRB wind to energy 
> 10^° eV, provided that T > 100 and that the magnetic field is close to equipartition 
with electrons. The former condition, F > 100, is remarkably similar to that inferred based 
on 7-ray spectra. There is no theory at present that allows a basic principles calculation of 
the strength of the magnetic field. However, magnetic field close to equipartition, ~ 1, 
is required in order to account for the observed 7-ray emission. 

We have assumed in the discussion so far that the fireball is spherically symmetric. 
However, since a jet-like fireball behaves as if it were a conical section of a spherical 
fireball as long as the jet opening angle is larger than T~^, our results apply also for a 
jet-like fireball (we are interested only in processes that occur when the wind is ultra- 
relativistic, F ~ 300, prior to significant fireball deceleration). For a jet-like wind, L in 
our equations should be understood as the luminosity the fireball would have carried had 
it been spherically symmetric. 

3. Implications of afterglow observations 

In addition to providing support to the validity of the qualitative fireball scenario 
described in §2.1, afterglow observations provide quantitative constraints on fireball model 
parameters. The determination of GRB redshifts implies that the characteristic GRB 7- 
ray luminosity and emitted energy are ~ lO^^erg/s and ~ lO^^erg respectively, an 
order of magnitude higher than the values assumed prior to afterglow detection. Afterglow 
observations also indicate that ~ ~ 0.1. This suggests that the constraint (|I|) is 
indeed satisfied, allowing proton acceleration to > 10^° eV. 

The observed GRB redshift distribution implies a GRB rate of -Rgrb ~ 10/Gpc^yr 
at z ~ 1. The present, z = 0, rate is less well constrained, since most observed GRBs 
originate at redshifts 1 < z < 2.5 Present data are consistent [§| with both no 
evolution of GRB rate with redshift, and with strong evolution (following, e.g., the lu- 
minosity density evolution of QSOs or the evolution of star formation rate), in which 
-Rgrb(^ = 1)/-Rgrb(^ = 0) ~ 8. The energy observed in 7-rays reflect the fireball energy 
in accelerated electrons. If shock accelerated protons and electrons carry similar energy, 
as indicated by afterglow observations, then the z = rate of cosmic-ray production by 
GRBs is similar to the generation rate of 7-ray energy, 

E\dhcR/dE),=o ^ lO^^erg/MpcVr, (3) 

where ( is in the range of ~ 1 to ~ 8. 

4. Comparison with UHECR observations 



In Fig. 1 we compare the UHECR spectrum, reported by the Fly's Eye the Yakutsk 
]T0| , and the AGASA |jll| experiments, with that predicted by the GRB model. The 
flattening of the cosmic-ray spectrum at ~ 10^^ eV, combined with the lack of anisotropy 
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Figure 1. The UHECR flux expected in a cosmological model, where high-energy protons 
are produced at a rate {E^dhcR/dE)z=o = 0.8 x lO^'^erg/Mpc^yr as predicted in the 
GRB model [Eq. (|)] (The flux above 10^^ eV is not sensitive to the z dependence 
of E'^dricR/dE). la flux error bars are shown. The highest energy points are derived 
assuming the detected events (1 for Fly's Eye and Yakutsk, 4 for AG AS A) represent a 
uniform flux over the energy range 10^° eV-3 x 10^° eV. 



and the evidence for a change in composition from heavy nuclei at low energy to light 
nuclei (protons) at high energy |]12|, suggest that an extra-Galactic source of protons 
dominates the flux a.t E > 10^^ eV. The UHECR flux predicted by the GRB model is in 
remarkable agreement with the observed extra-Galactic flux. 

The suppression of model flux above 10^^ '' eV is due to energy loss of high energy 



protons in interaction with the microwave background, i.e. to the "GZK cutoff" [|13 
Both Fly's Eye and Yakutsk data show a deflcit in the number of events, consistent 
with the predicted suppression. The deflcit is, however, only at a 2a confldence level 
[^]. The AGASA data is consistent with Fly's Eye and Yakutsk results below 10^° eV. 
A discrepancy may be emerging at higher energy, > 10^° eV, where the Fly's Eye and 
Yakutsk experiments detect 1 event each, and the AGASA experiment detects 6 events 
for similar exposure. 

The flux above lO^^eV is dominated by sources at distances < 40 Mpc. Since the dis- 
tribution of known astrophysical systems (e.g. galaxies, clusters of galaxies) is inhomo- 
geneous on scales of tens of Mpc, signiflcant deviations from model predictions presented 
in Fig. 1 for a uniform source distribution are expected above 10^°eV. Clustering of 
cosmic-ray sources leads to a standard deviation, a, in the expected number, N, of 
events above 10^° eV, given by a/N = 0.9((io/10Mpc)°'^, where do is the unknown scale 
length of the source correlation function and d^ ~ 10 Mpc for fleld galaxies. 

An order of magnitude increase in the exposure of UHECR experiments, compared to 
that available at present, is required to test for the existence of the GZK cutoff. Such 
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exposure would allow this test through an accurate determination of the spectrum in 
the energy range of 10^^''^ eV to 10^° eV, where the effects of source inhomogeneities are 



expected to be small |jT^. Moreover, an order of magnitude increase in exposure will also 
allow to determine the source correlation length do, through the detection of anisotropics 
in the arrival directions of ~ 10^^'^ eV cosmic-rays over angular scales of ~ c/q/SO Mpc 



5. GRB model predictions for planned UHECR experiments 

The rate at which GRBs occur within a distance of ~ 100 Mpc from Earth, the distance 
to which > 10^° eV proton propagation is limited due to interaction with the microwave 
background, is ~ 1 per 100 yr. This rate can be reconciled with the detection of several 
> lO^'' eV events over a period of a few years only if there is a large dispersion, > lOOyr, in 
the arrival time of protons produced in a single burst. The required dispersion is likely to 
result from deflection by random magnetic fields 0]. A proton of energy E propagating 
over a distance D through a magnetic field of strength B and correlation length A is 
deflected by an ang le 9, ~ {D/XY'^X/Rl, which results in a time delay, compared to 
propagation along a straight line, t{E,D) ^ 6'^D/Ac oc B^X. The random energy loss 
suffered by > 10^°eV protons as they propagate, owing to the production of pions, implies 
that protons observed at Earth with given energy have different energy histories along 
their propagation path. Thus, magnetic field deflection results not only in a delay, but 
also in a spread in arrival time of protons of fixed energy, comparable to the delay r. 

The current upper bound on the inter-galactic mag netic field 0, EA^/^ < iq-9q Mpc^/^ 
allows a spread t{E = 10^°eV, D = lOOMpc) ~ 10^ yr, well above the minimum, 
r ~ 100 yr, required in the model. The magnetic field upper bound implies an up- 
per bound on the number of GRBs contributing to the > 10^° eV flux at any given time, 
~ 10^/100 = 10^. The upper bound on the number of sources contributing to the flux 
above E decreases rapidly as E increases beyond 10^° eV, as the propagation distance 
decreases with energy due to the increase in pion production energy loss rate. This rapid 
decrease implies that at ii^ ~ 3"^^ eV there can be only a few sources contributing to the 



flux at any given time 

The GRB model therefore makes a unique prediction |16]: The UHECR flux at energy 
ii^ > 3 X 10^° eV should be dominated by a few sources on the sky. These sources should 
have narrowly peaked energy spectra, and the brightest sources should be different at 
different energies. This is due to the fact that at any fixed time a given burst is observed 
in UHECRs only over a narrow range of energy: If a burst is currently observed at some 
energy E then UHECRs of much lower (higher) energy from this burst will arrive (have 
arrived) mainly in the future (past). Testing the GRB model predictions requires an 
exposure 10 times larger than that of present experiments. Such increase is expected to 
be provided by the planned Auger |T^ detectors. 



6. 10 eV Neutrinos 

Protons accelerated in the fireball to high energy lose energy through photo-pion pro- 
duction in interaction with fireball photons. The decay of charged pions results in the 
production of high energy neutrinos [|l8l. The observed energy of a proton, for which the 
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observed 1 MeV photons are at the threshold of the A-resonance, is 0.2GeV^r^/l MeV. 
Typically, the neutrino receives ~ 5% of the proton energy. Thus, the typical energy of 
neutrinos resulting from interaction of accelerated protons with GRB photons is 

Et^5x 10''Tl,,eV. (4) 

The flux normalization is determined by the efficiency of pion production. The fraction 
of energy lost to pion production by protons producing the neutrino flux above is 



essentially independent of energy and is given by [jT8 

/. - 0.2-^^. (5) 

300^''10ms 

If GRBs are the sources of UHECRs, then using Eq. (|^) and and Eq. (j^) with C — the 



expected GRB neutrino flux is |18 




El^,^ ^ 1.5 X 10'^ {t;^]^ min{l, E,/Et}-^-, (6) 



GeV 
cm^ssr 

where stands for i/^, i/^ and Ug. 

The flux of ~ 10^^ eV neutrinos given in Eq. (P) implies that large area, ~ Ikm^, high- 
energy neutrino telescopes, which are being constructed to detect cosmologically distant 
neutrino sources |]T9| , would observe several tens of events per year correlated in time and 



in arrival direction with GRBs. Detection of neutrinos from GRBs could be used to test 
the simultaneity of neutrino and photon arrival to an accuracy of ~ 1 s (~ 1 ms for short 
bursts), checking the assumption of special relativity that photons and neutrinos have the 
same limiting speed to one part in 10^^, and the weak equivalence principle, according to 
which photons and neutrinos should suffer the same time delay as they pass through a 
gravitational potential, to one part in 10^ (considering the Galactic potential alone). 

The model discussed above predicts the production of high energy muon and electron 
neutrinos. However, if the atmospheric neutrino anomaly has the explanation it is usually 
given, oscillation to i/^'s with mass ~ 0.1 eV ||2D|, then one should detect equal numbers 
of j/^'s and z/^'s. Up-going r's, rather than /i's, would be a distinctive signature of such 
oscillations. Since z/^'s are not expected to be produced in the flreball, looking for r's 
would be an "appearance experiment." To allow flavor change, the difference in squared 
neutrino masses, Am^, should exceed a minimum value proportional to the ratio of source 
distance and neutrino energy. A burst at 100 Mpc producing lO^^eV neutrinos can test 
for Am^ > 10~^^eV^, 5 orders of magnitude more sensitive than solar neutrinos. 
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